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ABSTRACT 


^A spectroscopic investigetlon of a plasma shock tube may reveal 
the number density and temperature of the constituents of the plasma. 

Such an Investigation was conducted on a cylindrical fused quartz shock 
tube and a square, three Inch aluminum shock tube st thi rhyslsal— 
ta borato f »F- Aa t a sp aca-Corporati on. h oa -Anaftlaav -Caltfernla In June-July, 
1961. Deuterium and hydrogen plasmas were Investigated separately. Gen¬ 
eral characteristics Including ringing frequency and longitudinal visible 
shock front velocity were determined. Spectral line analysis from a 1.7 
meter quartz prism spectrograph permitted the determination of the shock 
tube Impurities and relative line Intensities. Use of the latter instru¬ 
ment in conjunction with a tungsten ribbon calibration lamp and a step 
filter provided an absolute intensity standard. For the hydrogen plasma, 
relative and absolute spectral line intensities werg,obtained whlch^er- 
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1. Introduction 

In recent years, particular attention has been focused on spectro¬ 
scopic diagnostics of plasma devices. Early Investigations and prelim¬ 
inary theories were based on the results of flame, arc, and spark ex¬ 
periments. Refinement to these theories and formulation of new theories 
have emerged from sophistication of the earlier plasma devices to Include 
among others that of the plasma shock tube. Within such a shock tube, 
the plasma Is formed by transferring energy from a shock generator to 
the gas filling. The presence of the shock wave and subsequent energy 
transfer Is evident from the appearance of a luminous shock wave within 
the tube, the properties of this wave being the subject of the spectro¬ 
scopic studies. 

The two spectroscopic Instruments used In this shock tube Investi¬ 
gation were the streak camera and the spectrograph. From these Instru¬ 
ments, the derived experimental data permitted determination of the 
following plasma parameters: 

a. shock tube Impurities 

b. luminous shock front velocity 

c. electron and Ionic temperatures 

d. electron and ion number densities 

Specific theories used to calculate these parameters Include those of 
relative and absolute line Intensities, line broadening, and apparent 
series limit. Throughout this paper, we have briefly mentioned each 
theory prior to applying the pertinent results derived. A more complete 
explanation of these theories will be found In the references provided. 

The actual Investigation was conducted at the Aerospace Corporation 
Physical Research Laboratory on two separate shock tubes, each briefly 
described In the next section. Calculations were Initiated at Aerospace 
Corporation and completed at the U.S. Naval Postgraduate School. 

It should be pointed out that not all of the previously mentioned 
plasma parameters were determined for each shock tube. Our study Is 
only Intended as a preliminary Investigation to future additional and 
more extensive studies to be performed on both shock tubi s. 
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2. Experimental Set-up 

Components comprising the experimental set-up may be generally div¬ 
ided into shock tubes, magnetic piston driver and related power components, 
gas supply and vacuum system, and monitor and test equipment. A photo¬ 
graph, Fig. 1, shows the physical placement o£ the experimental set-up 
during a spectroscopic investigation. A schematic of the experimental 
setup is shown in Fig. 2. In this section, only a brief description will 
be given of the components mentioned above. 

A. Shock tubes 

The first shock tube placed under investigation was a General 
Electric, type 204, fused quartz tube. Approximate dimensions 
were 45 inches in length with Inside and outside diameters of 
74 and 79 mm, respectively. A chemical analysis by the manu¬ 
facturer was conducted to determine tube material impurities 
(see Table 1). Prior to making a spectroscopic investigation, 
the tube was modified by the placement of two replaceable 
one-inch diameter flat optical quality sapphire windows about 
10 and 20 centimeters from the driver assembly. The latter 
modification prevented the necessity of making spectral ob¬ 
servations through cloudy tube inner surface impurity film 
deposits. 

The second shock tube was rectangular in shape and was fabri¬ 
cated of aluminum by Aerospace Corporation. Approximate di¬ 
mensions were an overall length of 36 inches and an inner 
and outer dimension of three and four inches square respective¬ 
ly. Spectroscopic observations were made through six inch 
diameter optical quality fused quartz windows 12.6 inches from 
the driver assembly. The smaller windows as noted in Fig. 3 
closer to the driver assembly were for later investigations 
(Interferometric) to be conducted on the shock tube. 

B. Magnetic Piston Driver and Related .'ower Components 

The magnetic piston driver was a modified form of the "conical 
shock tube" type^. Power for the magnetic field driver was 
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derived from a 40 KV potential between a capacitor bank of 
eight General Electric (type 14F97) 1.0^ £ can capacitors 
(25 KV) and a floating ground strap at a negative 15 KV 
potential (see Fig. 2.). 

C. Gas Supply and Vacuum System 

The vacuum system consisted of a fore pump, diffusion pump, 
and liquid nitrogen trap, an arrangement providing pressures 
as low as 10 ^ cm Hg. Vacuum gauges used were a Veeco, type 
RG-2A (pressures less than 100>4Hg) and a Autovac, type 3294B 
(pressures above lOO^Hg). In addition, a purge pump was in¬ 
corporated into the gas supply line. The gas supply system 
consisted of bottled gas and associated tubing and valve com¬ 
ponents. 

D. Monitor and Test Equipment 

The major components of moniter and test equipment were: 

1) Gaertner Scientific Corporation Quartz Spectrograph, 
type L 254 

2) Philips Research Laboratories Tungsten Ribbon Lamp, 
type U2 (with intensity calibration) 

Additional equipment used in data reduction included: 

1) Leeds and Northrup Recording Microphotometer, type 
6/00-2 

2) McPherson Instrument Measuring Comparator, type 100 
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3. Invc8t<g»tion Procedures 

A brief discussion Is made below outlining the preliminary tests made 
cn the shock tube system. Mention is then made of the procedures followed 
In the subsequent spectroscopic investigation. 

A. Preliminary 

The ringing frequency of the discharge circuit was determined 
by placing a pickup coil about one-half inch above the return 
straps of the driver assembly. Analysis of an oscilloscope 
Polaroid picture of the coll signal revealed a discharge 
circuit ringing frequency of 245 KC. 

The luminous shock wave produced in the fused quartz shock 
tube was examined with a rotating mirror streak camera. 

Fig. 4 is an enlargement of a typical luminous shock wave 
velocity smear picture. The appearance of weak shock waves 
after the primary brlf "it shock may be attributed to the 
ringing of the discharge circuit after firing. The heavy 
dark lines in Fig. 4 are reference distance markers along 
the shock tube axis. The entrance slit of the streak camera 
was aligned along the central shock tube axis in order that 
the 35 mm film strip yielded a plot of light intensity distri¬ 
bution along the tube axis versus time. 

The velocity of the luminous shock front was determined as a 

function of initial Hydrogen and Deuterium inlet pressure so 

!• 

that the optimum inlet pressure for maximum shock velocity 

* 

could be obtained . Fig. 5 and 6 are curves of initial in¬ 
let pressure versus shock velocity for Hydrogen and Deuterium 
gases respectively. The optimum pressures for maximum shock 
velocities were determined to be 8(^ of Hydrogen and of 

Deuterium yielding maximum shock velocities of 17.22 cm/^ sec 
and 11.65 cm/yOt^sec respectively. 

* The method used to calculate the velocity of the luminous shock front 
is shown in Appendix 1. 
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B. Spectroscopic 

Prior to firing either shock tube, the system was evacuated 
to about 10 cm Hg pressure. The gas supply line was then 
purged, after which the gas to be ionized, hydrogen or 
deuterium, was flushed through the system at approximately 
100 microns pressure for periods in excess of 15 minutes. 

The gas pressure was then regulated to that pressure corree- 
ponding to maximum luminous shock front velocity and the tube 
fired. Between tube firings, flushing occurred during which 
time the system returned to the desired pressure without any 
readjustment of the gas supply system. 
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4 Relative and Absolute Line Intensities 

Both relative and absolute line intensities are utilized in the 

spectroscopic calculation of electron temperature and number density. 

A tungsten ribbon standard lamp exposure on the same spectroscopic plate 

with the experimental exposures provides the reference for absolute and 

relative intensity calculations. Densitometer traces were made of the 

standard lamp and experimental exposures. Accepted data reduction pro- 
2 3 4 

cedures ’ * were then used in the calculation of relative and absolute 
line intensities. 

A. Relative Line Intensity 

In the calculation of relative intensities, the spectro¬ 
scopic plate was calibrated using the standard lamp source 
in order to account for the change in emulsion sensitivity 
with wavelength and the non-linear dispersion of the spectro¬ 
graph used. A calibatlon curve was then plotted, Fig. 7, 
which relates observed line intensities at all wavelengths 
to a common intensity scale as follows: 



BCAP 1 

h 1 

observed 

^2 

corrected ® ^^2^ 1 

^a 1 


(4-1) 


vAiere B(A.) is the correction factor which depends on 
emulsion sensitivity and spectrograph dispersion. 
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The 

corrected relative 

line Intensities 

were: 

Element 

Multiplet^'^ 

Wavelength 

Relative Intensity 




normalized to 



6563 

1.00 



A861 

0.0282 

Hr 


A3A0 

0.00905 



A102 

0.00370 



3970 

0.00155 

0 II 

3 

37A9 

0.0211 

0 II 

36 

A190 

0.003A6 

0 II 

2 

A3A9 

0.0106 

0 II 

5 

AA15 

0.00780 

0 II 

1 

A6A9 

0.0150 

Ca II 

A 

3179 

0.0A55 

Ca II 

3 

3737 

0.00308 

Ca II 

1 

393A 

0.08A3 

Ca II 

1 

3968 

0.0575 

Mg II 

1 

2796 

1.02 

Mg II 

1 

2803 

0.906 

Mg II 

3 

2798 

0.262 

Mg II 

2 

2936 

0.12A 

Ti II 

2 

3235 

0.105 

T1 II 

2 

32A2 

0.0555 

Tl II 

1 

33A9 

0.118 

Na I 

6 

5688 

0.00628 

Na I 

1 

5890 

0.360 

Na I 

1 

5896 

0.310 

Na I 

5 

6161 

O.OIAI 1 

C II 

6 

3921 

O.OOAIO 

C II 

A 

A267 

0.037A 1 
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B. Absolute Line Intensities 

The absolute line intensities for the Balmer series lines 
were calculated separately for through H£- using cali¬ 
brated standard lamp exposures for each of these wavelengths 
in the general method described by Harrison, Lord, and 
Loofbourow^. 

Several assumptions were made: 

1) In order to obtain a value for the spectroscopic 

system transmission factor, the following representa- 

4 

tive transmission values were assumed : 

a. 0.90 for each lens 

b. 0.85 for each prism 

c. 1.00 for each mirror 

In order to minimize the error Involved in this assumption, 
the standard source was used with the same optical parameters 
that were used in the experiment^. 

2) The path length from the standard lamp to the entrance 
slit was sufficiently short to avoid appreciable hori- 

7 

zontal extinction of the standard source radiant energy . 

3) To obtain the radiant energy, the solid angle subtended 

4 

by the collimator at the slit was approximated by : 

W = ^ (4-2) 

F 

where W = solid angle subtended by the collimator 
at the slit 

A = area of the collimator 
F = collimator focal length 

4) The vertical variation of the emulsion was negligible. 

5) The time of radiation was approximately 5 x 10 ^ seconds. 

This time corresponds to the duration of the bright 

region in the smear photograph. Fig. 3. 

0 

6) The shock tube was optically thin . 
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Tht absolute energies and intensities calculated per shock 
tube firing were: 


Line 


“6 


Energy 

Intensity 

(ergs) 

, -2 -1 , -1. 
(ergs-cn -sec -ster ) 

1.353 X lO"^ 

0.6950 X 10^ 

0.386 

0.1980 

0.133 

0.0683 

0.0288 

0.01478 

0.0220 

0.01130 
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5. Shock Tube Impurities 

The Identification of shock tube Impurities for both shock tubes, 
although a tedious procedure, was of primary Interest and produced a 
useful result. Following the extension of the spectrograph calibration 
curve to Include the 2000-7000 X range, shock tube firings were conducted 
under the conditions shown In Table II. In conjunction with the above 
swntloned curve and comparator measurements, line wavelengths were obtained. 
A typical 3400 center wavelength spectrum Is shown In Fig. 8. Element 
and sniltlplet Identification was then made using tables prepared by 
Charlotte Moore^* . Line Identification results are shown In Table III 
and are valid for both the fused quartz shock tube and fused-quartz-wlndow 
aluminum shock tube. This Impurity spectra corresponds to the known con¬ 
stituents of the shock tube electrodes and walls and verifies that ablation 
occurs during shock tube firing. This result will prove useful In Impurity 
Identification In similar shock tube systems employing fused quartz sur¬ 
faces as spectral windows. 
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6. Electron Number Densities 

A knowledge of number densities is necessary in several of the methods 
for determining electron temperature. Two approaches to number density 
calculation are those of A. apparent series limit and B. the Stark broadening 
profile comparison. The formulation of a hydrogen plasma clearly validates 
the assumption of equating electron and ionic number densities. 

A. Apparent Series Limit 

As the principal quantimi number, m, of a Balmer line increases. 
Stark broadening of that line increases. Therefore, these 

9 

lines merge and form an apparent continum. Inglis and Teller 
derived a relationship between the last discernible line before 
the resulting apparent continum and number density for ap¬ 
plication to astrophysical observations. A slight modification 

Q 

to their equation was made by Lochte-Holtgreven for use with 
radiating plasma observations. This relationship is: 

Log N ■ 23.46 - 7.5 Log m (6-1) 

max 


where N a n = n.anumber density 
6 x> 

principal quantum number of the last 
discernible line within the series. 


m 


max 


Letting ^ for IL (see Table Ill)re8ult8 in the valuee: 

1 -3 

n 3 n, = 1.3 X 10 cm 
e i 


6. Stark Broadening Profile 

For particularly broad Balmer line prof ilea, it may be aaaum'ed 
linear Stark broadening is dominant. Among the Initial treat¬ 
ments of this phenomenon was that of Holtzmark, who considered 

o 

only the effects of perturbing ions on the emitting atoms . 

As this theory did not closely agree with experimental line 
profile observations, further refinement of this theory was 
formulated by Griem, Kolb, and Shen and included the per¬ 
turbing effects of both electrons and ions on the hydrogen 

and are shown 


atoms. The experimental line profiles of 
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In Fig. 9 and 10, each a plot of normalized line intensity, 

I/I , versus displacement from center wavelength,A A* 

° 10 
Theoretical profiles, as determined by Griem, Kolb, and Shen 

are also shown superimposed on the experimental profiles. In 

plotting the theoretical profiles, the useful relationships 

developed were : 

AA= 1-25 X 10® a (6-2) 


I/I = sCg) 


0.625 X lO"® 


(6-3) 


where N > electron number density. 

a s ratio of the distance from the unperturbed 

line, A A, to the Holtzmark normal field 

.2/3 

S(a) s line profile function such that |s(C2) da ^ 1. 

-c-o 


strength, = 2.61 eN 


• oo 


(area) » area under the normalized experimental line 
profile. 


Reference to Fig. 9 and 10 indicates that the effect of tern- ' 
perature is small and that of number density is predominant in 
determining the theoretical profile. From the curve fitting 
then, it can be seen that 

10^^ < N < 10^^ 

where N = n = n, in centimeters 
e 1 

In the calculations of electron temperature to follow, it will 

be assumed that electron and ionic number densities equal 10^^ 
“3 

centimeters . This approximation is borne out by the above 
calculations and is considered accurate within an order of 
magnitude. 
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7. Electron Temperature Determination 

The temperature of a pLisma may be defined only if there is statistical 
11 12 

or thermal equilibrium * . Assuming thermal equilibrium, the population 

distribution over the excited energy levels varies according to Boltzmann's 
law. The spectral intensity may therefore be related to the temperature. 

The intensity of an emission line corresponding to a transition from 
n to m is given by: 


I 

nm 


B W™' hU 


nm 


BA™" H hi/ 

n nm 


(7-1) 


where B = instrumental constant 

W™" = N A™ = number of transitions taking place in light source 

second 


N 


.nm 


hU 


nm 


= number of atoms in initial state n 


= Einstein transition probability = 
= radiated photon energy 


number of transitions of atom 
sAcond 


For the case of thermal or statistical equilibrium, the population 
distribution may be defined as: 

"a " -T «'■ 


where N « total density of emitting species in volume for which the 
intensity is being measured, 
g^ = degeneracy of initial state n (statistical weight) 
u = partition function of all statistical weights 
= energy of the initial state n above ground state 
T^ = election temperature 

Also^^, g^ A™"= 2/3 X 10^^ (7-3) 

where J = 1 or 2 

wavelength of the transition 
f^a oscillator strength determined for absorption 

f 2 = oscillator strength determined for emission 
g^=g 2 =statlstical weight of the Initial state n 
gj^a statistical weight of the final state m 
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(7-4) 


Therefore, 

where D = (2^1 (2/3 x lo‘‘) ) 

c 

Equation 7"4 may be used with the absolute Intensity of a single line 
in order to determine the electron temperature, but this method is diffi¬ 
cult due to the uncertainty in the total density, M, of the plasma. In 

calculating N, it is necessary to determine the density of each constituent 

a 

present, requiring additional assumptions Lo those already madt . Also, 
there is uncertainty of reproducing the same constituent densities with 
each tube firing due to the vari'ible effects of ablation in the shock tube.^ 


Several methods were attempted in the determination of electron 
temperature. 

i-iethod I: Relative Intensity 

Equation 7-4 may be rearranged as follows: 

(7-5) 


Log 




E^ log e 
n 

kT 


Since D is independent of the energy E , it may take any value and 

n 

the left hand side of 7-5 may be expressed in arbitrary units . in the 

calculation of temperature using the lines from an element, the relative 

13 14 

intensities and relative values of (g f ) may be used ' . A plot of 


log 



versus E for a number of lines from an element should therefore result 

“ 1 
in a straight line with slope = - Log e. 
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Method II: Relative Intenaity Ratios 

Equation 7-4 may also be formed Into a ratio of two lines from the 
same element, giving the resultant equation: 



(7-6) 


Since ratios are being used any atbltrary consistent units may be 
used for I,y and (gf). A plot of 


log 




^ 2 ^ * 1^1 


versus (C^ - E^) for a number of pairs of lines from an element should 
therefore result in a straight line with slope = - log e . 


In Methods I and II, the resultant graphs, a sample shown in Fig. 11, 
not only yield a value for the electron temperature, but also ac' as a 
quasi-check for the assumptions of thermal equilibrium. If the resultant 
points do not fall on s. straight line and the data has been checked for 
systematic errors, then any curvature indicace» a lack of thermal equili¬ 
brium in the population of the excited states of the constituents and no 

11 13 

temperature can be defined (this curving was not observed in our results). * 


Method III Absolute Intensity 

The absolute intensity of the Balmer lines in emission may be used in 
the electron temperature determinatlcu as follows:^^ 


abs 




exp 


-E -f-E \ 
n o 

kT 


(7-7) 
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where e « electron charge 
h = Flunk's constant 
£ s oscillator strength 

n s principle quantum number of the upper state 

= neutral atom number density 

1 3 radiating layer thickness 

m 3 electron mass 
e 

3 emission line wave length 

s energy of the state n 

= energy of the ground state = - 13.595 ev 


The depression of ionization potential within the plasma is a result 
of the presence of ions which are a finite distance apart. The effect 
of the ions is to alter the Coulombic field in which the electrons move 
and may be expressed as:^ 


AE = 7.0 X 10“^ ev (7-8) 

where = electron number density. 

17 -3 

This corresponds to a depression of 0.32 ev with » 10 cm . 

Assuming that the thickness of the radiating layer is approximately 
given by the shock tube dlameter^^ and substituting the electron number 
density determined earlier, equations 7-7 and 7-8 reduce to: 

-1.16 X 10^(13.28 + E^) 

36.77 * A/l.^, am 

1 ‘ J 

where T is in degrees Kelvin, N in Centimeters) Ain centimeters, I . 

-2 “1 -1 ° 

in ergs-cm -sec -ster ., E^ in electron-volts, and f and n are dimen¬ 
sionless. 

Equation 7-9 may be solved for any of the Balmer lir^b in emission if 
the neutral atom number density can be obtained. This latter quantity was 
assumed equal to the electron number density, that is 50 percent ionization. 
The effect of this assumption is discussed in the next section. 
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Rcfults of the electron temperature calculations are 


Method 

Element 

Temperature 
degree Kelvin 

e-v 

I 

H 

6,300 

0.54 

I 

Na I 

50,000 

4.30 

I 

Ca II 

10,800 

0.93 

I 

C II 

31,300 

2.70 

1 

Mg II 

24,800 

2.14 

I 

Ti II 

6,600 

0.57 

II 

H 

6,800 

0.59 

II 

Na I 

47,000 

4.05 

II 

0 II 

24,000 

2.07 

III 

«a 

7,910 

0.68 

III 


7,890 

0.68 

III 

Hr 

7,730 

0.67 

III 

"(5 

7,300 

0.63 

III 


7,300 

0.63 
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8. Discussion 

The calculation of accurate temperatures, power levels, and densities 
in a plasma is very important. Aside from the ability to describe accurately 
the physical parameters of the plasma system, these values enable the further 
derivation of experimental transition probabilities and the calculation of 
collision cross-sections which to date are not sufficiently known. 

In association with Doctor Richard M. Head, considerable time was 
devoted to identification of the various impurities present in the shock 
tubes^^. Since the observed impurity spectrum was Identical for the quartz 
shock tube and for the aluminum shock tube with quartz windows, it is felt 
that the list of spectral lines shown in Table III will be of assistance in 
'' ipurity identification for most quartz systems operating under similar con- 
d'^tlons. It should be pointed out that only those elements that could con¬ 
tribute significantly to the intensity of a particular line are included in 
Table III. There were, in addition to the listed wavelengths, some 100 lines 
that were of very low intensity which are omitted due to uncertainty in their 
origin. Also, most of the lines below 2800 are omitted because of limited 
identification data in this region. 

Using the spectral lines of the various impurities in the calculation 
of temperature, a check may be made on whether a unique temperature can be 
defined for the plasma'*' . 

The accuracy of the calculated temperature and number densities in the 
shock tube is highly dependent on the relative and absolute line intensities 
obtained. In order to determine correctly the relative or absolute intensity 
of a line, a number of factors must be carefully considered. These Include: 

1) Spectroscopic techniques used 

2) background continimi correction 

3) plate-calibration corrections for 

a) dispersion 

b) eiulslon sensitivity 

c) standard source 

4) self-absorption considerations 

5) computational errors 
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The first three correction factors uere determined using standard 
spectroscopic techniques described In Harrison, Lord, and Loofbourow , 
and with the special techniques mentioned In appendix II. 

The problem of self-absorption. If present, can seriously affect the 
measured Intensities of the lines, and therefore lead to meaningless re¬ 
sults. Severe self-absorption Is characterized by flat-topped lines or 

separation of one line Into two lines by strong absorption In the middle 

13 

of the spectral line . For plasma systems operating at temperatures 

g 

greater than approximately one ev, the problem of self-absorption is lessened; 

Such absorption Is related to the total Intensity of emission from an element, 

which Is temperature (energy) dependent, and the total number of potential 

absorbers. At temperatures (energies) near the maximum Intensity of emission, 

the amount of absorption becomes negligible since the number of absorbers 

present Is much less than the number required to absorb all of the emitted 

energy. For hydrogen, this temperature corresponds to approximately one ev. 

As severe self-absorption was not observed, the system was therefore assumed 

to be optically thin for the Balmer series. Future experimental work with 

this system should Include a check on the validity of the optically thin 

plasma assumption by placing a mirror on the optical axis behind the source, 

13. 

as recommended by Pearce 

Since neither radial distribution not time resolution studies were con¬ 
ducted In this experiisent, the intensities determined permit only the cal- 

14 

culatlon of a mean integrated temperature for the plasma . 

Taking into account the above considerations. It Is believed that the 
intensities calculated are representative of the populations of the various 
energy states and consequently may be used In temperature and number density 
calculations using the Boltzmann distribution. 

The two methods, apparent series limit and Stark Broadening profile 
comparison, employed to determine number densities are based on the exist¬ 
ence of line broadening phenomena. They are, however, unrelated as to the 
paramet <rs being considered. 

In the method of the apparent series limit, there Is no known limita¬ 
tion placed on equation 6-1 other than a decrease In the calculated value 
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by a factor of two in number density If the plasma temperature is less 

5 8 

than 10 / m (where m = 7 in our study) . This seems unlikely from 
max max 

the temperature calculations in this research. 

In the method of Stark broadening profile comparison, the range of 

14 “3 19 

applicability for experimental profiles provided is 10 Z. n cm Z. 10 

— e — 

and 10,000<T ^40,000. As previously mentioned, the effect of number 

density variation is dominant over that of temperature variation. The 

accuracy of the theoretical profiles is approximately ten percent ^ iind 

the accuracy of experimental profiles used are estimated to be of this 

same order of magnitude. The experimental accuracy is not quite as high 

for succeeding lines in the Balmer series because of masking of the wings 

17 -3 

by impurities. The choice of 10 cm must be within an order of magnitude 

for electron number density since it appears consistent with the results of 
17 “3 16 “3 17 

1.3 X 10 cm and 10 <n^ cm <10 by the series limit and Stark broaden¬ 

ing methods respectively. 

Prior to electron temperature calculation using absolute intensities, 

it was necessary to estimate the neutral atom number density. The value 
17 -3 

10 cm was decided upon with the stipulation of studying its variation 
as a parameter. Variation of neutral atom number density by two orders 
of magnitude affected the resulting electron temperature by lees than a 
factor of two. In view of this result, it is not bellved that the elec¬ 
tron temperature results are materially affected by the necessity of this 
procedure. 

In the determination of electron temperature, the calculations were 

limited by the lack of transition probabilities for many of the observed 
18 

lines . In methods I and II the accuracy of the results depends strongly 
on the properties of the lines used in the calculations. It ia very de¬ 
sirable to choose a number of lines originating from a wide range of upper 

14 

energy levels and occurring over a small wavelength interval . This pro¬ 
cedure was somewhat compromised due to the lack of transition probabilities. 
The range of validity of methods I and II is up to approximately 5 x 10 \ 

g 

with the proper selection of the elements . For each wavelength line there 
is a relatively flat maximum in the temperature-intensity profile which ex¬ 
tends over a few thousand degrees. Therefore, when measuring a temperature. 
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It Is more accurate to select a line where the slope of the temperature- 
intensity profile is rapidly changing. From data presented by Lochte- 

g 

Holtgreven , the temperatures determined using Oxygen II and the Balmer 
lines were on the rapidly rising slopes of the temperature-intensity pro¬ 
files and thus in the more accurate range. Similar data were not available 
for the other elements. 

17 -3 

Since the electron density is relatively high (10 cm ), there is a 

rapid approach to local thermodynamic equilibrium. The formula for re- 

19 

laxation time is given by Spitzer : 


,~ 11.A ._3/2 

t . (T ) sec. 

N Z InA 


V 


where InA — 20 

As 1 for a hydrogen plasma 
Z s 1 for a hydrogen plasma . 

N n 10 cm 
T s 6270° K (0.54 ev) 

Therefore t —10 seconds or less, and the concept of electron temperature 

has some meaning in the usual statistical sense. This is further verified 

by the results of calculations used in methods I and II. In these methods 

a marked departure from equilibrium is characterized by a curving, or trend, 

13 

in the graphic data . This curving, and hence departure from equilibrium, 
was not observed in our result (see Fig. 11). Several other advantages were 
gained by using methods 1 and II. A greater number of lines could be used 
with these methods than in the absolute intensity method; and therefore, any 
anomally caused by an Improper f-value, self-absorption, or intensity error 
(which would result in a wide variation from the other valuer) could be easily 
distinguished. In addition the more easily measured relative intensity may 
be used instead of absolute intensity, and it is not necessary to know ac¬ 
curately the composition of the plasma. 
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The temperatures calculated by methods 1 through III varied from 
0.54 to 4.2 ev. This result is to be expected, since each plasma con¬ 
stituent is likely to exhibit a different temperature during the time 
Integrated observation. Future investigations should Include a tempera¬ 
ture-time variation study within the shock tubes. 

The results of method III using the absolute intensity agree closely 
with the electron temperature calculated using the relative intensities of 
the Balmer lines. The difference between the methods, which is less than 
a factor of 25 percent, shows that the assumptions made in method III and 
those used in the absolute intensity calculations were reasonable. 

The ion temperature may be estimated using the results of the absolute 

Intensity measurements. This problem has been investigated theoretically 

20 

for a pure hydrogen plasma considering the following factors : 

1) Energy transfer from atoms to electrons caused by elastic 
collisions between atoms and electrons and ions and electrons. 

2) The energy dissipation by the electrons that is determined by: 

a) inelastic collisions of electrons with atoms and electrons. 

b) the excess of emission over absorption in free-free transi¬ 
tions. 

c) the excess of emission over absorption in free-bound and 
bound-free transitions. 

The difference between electron and kinetic temperatures under these con¬ 
ditions would be: 

-3 

T,-T <10 T 
1 e e 

10 -3 

when N > 10 cm 
e' 

and 5000°K<T <10^OK 
e 

From this relation the ion temperature is approximately equal to the 

electron temperature for the conditions existing in the shock tube (l.e., 

N = 10^^ cm ^ and T > 0.54 ev.). 
e e — ' 
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Conclusions 

A. Shock Tube Characteristics 

Spectroscopic diagnostic techniques reveal the following plasma 
parameters existing within the two shock tubes with hydrogen 
Introduced; 

17 -3 

1) An electron and Ionic number density of approximately 10 ,cm , 

2) A peak temperature of approximately 4 ev. 

3) A condition of local thermodynamic equilibrium. 

B. General 

1) Shock tube systems employing quartz windows and operating under 
similar conditions are characterized by nearly Identical Im¬ 
purity spectra. 

2) Spectral line absolute Intensity measurements provide a useful 
mechanism for determining plasma temperature. 

3) Further Investigation of the shock tubes under consideration 
should Include time resolved and radial temperature and number 
density studies. 



Fig. 1. Experimental Set-up During Spectroscopic Investigation 
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Fig. 2. Schematic of the Experimental Set-up 





Fig. 3. Square Shock Tube 
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Fig. 4. Shock Wave Smear Photograph 




Wave Velocity, cm^Asec 



Fig. 5. Hydrogen Pressure Variation Effect on Shock Wave Velocity 
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Deuterium Pressure, microns 

F g. 6. Deuterium Pressure Variation Effect on Shock Wave Velocity 


29 



9*0 
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Fig. 7. Spectroscopic Plate Calibration 




Note: a-plate 1 of 7/18, b-plate 1 of 7/10; see Table II for firing conditions; exposures froa left 
to right - single shot, two shots superlnposed, -four shots superiaposed. 
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Impurity 

Parts per Million 

F.jOj 

12 

Ti Oj 

4 

^2 °3 

160 

Ca 0 

22 
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CM 
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Na^ 0 


Lij 0 

/ 9 

/ 
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0.94 


Table I • Fused Quartz Tube Impurities as 
Provided by Manufacturer 
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Table Ill Spectral Line Identification 


Wavelength 

(X) 

Element 

Multiplet^*^ 

Number 

Excitation Potential 
Low High 

Relative 

Intensity 

(visurA) 

2367.064 

A1 

I 

4 

0.01 

5.21 

10 

2395.06 

B 

II 

4 

12.64 

17.79 

6 

2435.160 

Si 

I 

45 

0.78 

5.85 

15 

2478.556 

C 

I 

61 

2.67 

7.65 

10 

21^96.773 

B 

I 

1 

0.00 

4.94 

18 

2497.724 

B 

I 

1 

0.00 

4.94 

20 

25l6.896 

Si 

I 

1 

0.01 

4.93 

18 

2514.315 

Si 

I 

1 

0.00 

4.91 

16 

2516.109 

Si 

I 

1 

0.03 

4.93 

20 

2519.203 

Si 

I 

1 

0.01 

4.91 

17 

2524.108 

Si 

I 

1 

0.01 

4.90 

18 

2528.510 

Si 

I 

1 

0.03 

4.91 

18 

2532.38 

Si 

I 

86 

1.90 

6.77 

8 

2541.83 

Si 

III 

8 

10.23 

15.09 

18 

2631.28 

Si 

1 

83 

1.90 

6.59 

12 

2790.768 

Mg 

II 

3 

4.40 

6.83 

4 

2795.523 

Mg 

II 

1 

0.00 

4.41 

22 

2797.989 

Mg 

II 

3 

4.41 

8.83 

9 

2802.698 

Mg 

II 

1 

0.00 

4.40 

21 

2816.189 

A1 

II 

7 

7.39 

11.77 

10 

2852.120 

Mg 

I 

1 

0.00 

4.33 

2 

2881.595 

Si 

I 

43 

0.78 

5.06 

24 

2904.29 

Si 

II 

17 

9.79 

14.04 

6 

2905.70 

Si 

II 

17 

9.80 

14.04 

10 

2936.96 

Mg 

II 

2 

4.41 

8.62 

1 

2987.65 

Si 

I 

1 

0.78 

4.91 

16 

3072.971 

Ti 

II 

5 

0.00 

4.02 

0 

3075.225 

Ti 

II 

5 

0.01 

4.02 

0 

3078.645 

Ti 

II 

5 

0.03 

4.04 

1 

3082.159 

A1 

I 

3 

0.00 

4.00 

12 

3086.225 

Si 

III 

1 

17.63 

21.63 

9 

3088.027 

Ti 

II 

5 

0.05 

4.05 

3 
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Wavelength 

Element 

Multiplet^’^ 

Bxcitat 

ion Potential 

Eelative 

(X) 



Number 

Low 

High 

Intensity 

(visual) 

3092.716 

A1 

I 

3 

0.01 

4.00 

17 

3092.8A3 

A1 

I 

3 

0.01 

4.00 

3096.786 

SI 

HI 

1 

17.64 

21.62 

1 

3129.76 

Zr 

II 

5 

0.04 

3.98 

0 

3129.16 

Zr 

II 

23 

0.52 

4.47 

L 

3129.44 

0 

II 

14 

25.53 

29.47 


3149.56 

SI 

IV 

2 

26.95 

30.86 

6 

3165.72 

Si 

IV 

2 

26.97 

30.86 


3165.98 

Zr 

II 

5 

0.16 

4.06 

12 

3166.29 

Zr 

II 

48 

0.80 

4.70 


3179.34 

B 

II 

Unci. 



1 

3179.332 

Ca 

II 

4 

3.14 

7.02 

4 

3203.89 

Si 

II 

7 

10.02 

13.87 

6 

3210.04 

SI 

II 

7 

10.03 

13.87 

12 

3210.52 

Si 

III 

Unci. 



3234.517 

Ti 

II 

2 

0.05 

3.86 


3234.00 

SI 

III 

6 

21.62 

25.44 

0 

3236.573 

Tl 

II 

2 

0.03 

3.84 

4 

3239.037 

Ti 

II 

2 

0.01 

3.82 

2 

3241.984 

Ti 

II 

2 

0.00 

3.81 


3241.67 

Si 

III 

6 

21.63 

25.44 


3247.540 

Cu 

I 

1 

0.00 

3.80 

20 

3261.596 

Ti 

II 

66 

1.23 

5.01 

0 

3261.596 

Ti 

II 

89 

1.88 

5.67 

L 

3265.46 

0 

III 

8 

36.32 

40.10 

1 

3273.957 

Cu 

I 

1 

0.00 

3.77 

19 

3279.26 

Zr 

II 

3 

0.09 

3.86 

4 

3282.329 

Ti 

II 

66 

1.22 

4.98 

1 

3284.72 

Zr 

II 

4 

0.00 

3.76 

1 

3302.34 

Na 

I 

2 

0.00 

3.74 


3302.94 

Na 

I 

2 

0.00 

3.74 

0 

3322.936 

Ti 

II 

7 

0.15 

3.86 

A 

3323.61 

B 

II 

Unci. 




3341.875 

Ti 

II 

16 

0.57 

4.26 

6 

3343.770 

Ti 

II 

7 

0.15 

3.84 

2 

3349.399 

Ti 

II 

1 

0.05 

3.73 

10 

3361.213 

Ti 

II 

1 

0.03 

3.70 

7 
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Wavelength 

Element 

Multiplet^’^ 

Excitation Potent ' il 

Bclative 

(X) 



Number 

Low 

High 

Intensity 

(visual) 

3372.800 

Ti 

11 

1 

0.01 

3.67 

5 

3383.761 

Ti 

II 

1 

0.00 

3.65 

4 

3388.29 

Zr 

II 

2 

0.00 

3.64 


3388.87 

Zr 

II 

74 

0.97 

4.61 


3391.96 

Zr 

II 

1 

0.16 

3.80 

9 

3410.26 

Zr 

II 

11 

0.41 

4.03 

1 

3430.53 

Zr 

II 

11 

0.46 

4.06 

1 

3438.23 

Zr 

II 

1 

0.09 

3.68 

8 

3451.41 

B 

II 

1 

9.06 

12.64 

24 

3452.28 

B 

II 

Unci. 



3479.39 

Zr 

II 

46 

0.71 

4.26 

2 

3481.14 

Zr 

11 

46 

0 

00 

• 

4.34 

3 

3483.54 

Zr 

II 

33 

0.75 

4.30 

1 

3496.18 

Zr 

II 

1 

0.04 

3.57 

7 

3505.67 

Zr 

II 

1 

0.16 

3.68 


3505.47 

Zr 

II 

90 

1.53 

5.05 

3 

3504.890 

Ti 

II 

88 

1.88 

5.40 


3510.840 

Ti 

II 

88 

1.88 

5.40 

1 

3535.40C 

Ti 

II 

98 

2.05 

5.54 

1 

3542.65 

Zr 

II 

113 

1.75 

5.24 

1 

3551.94 

Zr 

II 

1 

0.09 

3.57 

2 

3556.61 

Zr 

II 

9 

0.46 

3.93 

3 

3572.47 

Zr 

II 

1 

0.00 

3.45 

5 

3576.88 

Zr 

II 

9 

0.41 

3.86 

2 

3586.9 

A1 

II 

7(9)* 

11.80 

15.24 

15 

3590.46 

Si 

III 

7 

21.79 

25.22 

4 

3596.048 

Ti 

II 

15 

0.60 

4.04 

1 

3601.623 

A1 

III 

1 

14.31 

17.74 

5 

3612.352 

A1 

III 

1 

14.31 

17.73 

3 

3614.79 

Zr 

II 

9 

0.36 

3.77 

1 

3685.192 

Ti 

II 

14 

0.60 

3.95 

A 

3685.192 

Ti 

II 

14 

0.57 

3.92 


3698.17 

Zr 

II 

71 

1.01 

4.34 

2 

3713.103 

A1 

III 

4 

17.74 

21.07 

1 

* 9 lines contribute to this multiplet. 
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Wavelength Element Multiplet^’^ 
(a) Number 


3727.33 

0 

II 

3 

Zin.12 

Zr 

II 

112 

3736.901 

Ca 

II 

3 

3749.49 

0 

II 

3 

3751.60 

Zr 

II 

71 

3759.291 

Ti 

II 

13 

3759.87 

0 

III 

2 

3761.320 

Ti 

II 

13 

3791.41 

SI 

III 

5 

3791.26 

0 

III 

2 

3796.11 

Si 

III 

5 

3806.56 

Si 

III 

5 

3836.76 

Zr 

II 

16 

3853.657 

Si 

II 

1 

3856.021 

Si 

II 

1 

3862.592 

Si 

II 

) 

3881.97 

Zr 

II 

134 

3882.197 

0 

II 

12 

3905.527 

Si 

I 

3 

3911.960 

0 

II 

17 

3918.997 

C 

II 

4 

3920.677 

C 

II 

4 

3924.44 

Si 

III 

Unpl 

3933.664 

Ca 

II 

1 

3944.009 

A1 

I 

1 

3961.523 

A1 

I 

1 

3968.470 

Ca 

II 

1 

3970.074 



1 

3973.263 

0 

II 

6 

3991.14 

Zr 

II 

30 

3998.98 

Zr 

II 

16 

3998.635 

Ti 

I 

12 

4062.90 

0 

II 

50 

4069.897 

0 

II 

10 

4069.636 

0 

II 

10 


Excitation Potential 
Low High 

Relative 

Intenalty 

(visual) 

22.88 

1.74 

26.19 

5.05 

1 

3.14 

6.44 

I 

22.90 

26.19 

2 

0.97 

4.26 

1 

0.60 

33.04 

3.89 

36.32 

4 

0.57 

3.85 

3 

21.62 

33.04 

24.88 

36.29 

2 

21.62 

24.88 

4 

21.63 

24.88 

7 

0.56 

3.77 

2 

6.83 

10.03 

9 

6.83 

10.03 

24 

6.83 

10.02 

21 

2.40 

25.55 

5.58 

28.73 

1 

1.90 

5.06 

22 

25.55 

28.71 

1 

16.26 

19.41 

1 

16.26 

19.41 

1 



1 

0.00 

3.14 

18 

0.00 

3.13 

14 

0.01 

3.13 

16 

0.00 

3.11 

14 

10.15 

13.26 

IW 

23.34 

26.45 

1 

0.75 

3.85 

2 

0.56 

0.05 

3.64 

3.13 

1 

28.58 

31.62 

I 

25.53 

25.52 

28.56 

28.56 

3 
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Wavelength 

Element 

Multiplet 

(X) 

4072.164 

0 

II 

Number 

10 

4076.00 

C 

II 

36 

4075.868 

0 

II 

10 

4088.863 

SI 

IV 

1 

4097.260 

0 

II 

20,48 

4096.543 

0 

II 

21 

4101.737 

% 

Si 


1 

4102.926 

I 

2 

4116.104 

Si 

IV 

1 

4121.95 

B 

II 

2 

4128.053 

Si 

II 

3 

4130.884 

SI 

II 

3 

4149.22 

Zr 

II 

41 

4149.897 

A1 

III 

5 

4150.138 

A1 

III 

5 

4189.788 

0 

II 

36 

4267.27 

C 

II 

6 

4267.02 

C 

II 

6 

4275.52 

0 

II 

67 

4340.468 



1 

4349.426 

0 

II 

2 

4347.425 

0 

II 

16 

4366.896 

0 

II 

2 

4414.909 

0 

II 

5 

4416.975 

0 

II 

5 

4479.966 

A1 

III 

8 

4479.891 

A1 

III 

8 

4512.535 

A1 

III 

3 

4529.176 

A1 

III 

3 

4552.654 

SI 

III 

2 

4567.872 

Si 

III 

2 

4641.811 

. 0 

II 

1 

4649.139 

0 

II 

1 

4661.635 

0 

II 

1 

4861.332 



1 


Excitation Potential 
Low High 

Relative 

Inteniity 

(visual) 

25.54 

28.57 

3 

24.27 

25.55 

27.30 

28.58 

6 

23.95 

26.97 

8 

25.72 

25.73 

28.73 

28.74 

2 

10.15 

13.16 

3W 

1.90 

4.91 

2 

23.95 

26.95 

3 

18.60 

21.59 

5 

9.79 

12.78 

23 

9.80 

12.78 

24 

0.80 

3.77 

2 

20.47 

20.47 

23.44 

23.4A 

2 

28.24 

31.18 

2 

17.97 

17.97 

20.86 

20.86 

6 

28.73 

31.62 

i 

10.15 

13.00 

low 

22.90 

25.55 

25.74 

28.39 

2 

22.90 

25.73 

2 

23.34 

23.32 

26.14 

26.11 

3 

20.69 

20.69 

23.45 

23.45 

2 

17.73 

20.47 

1 

17.74 

20.47 

2 

18.92 

21.63 

4 

18.92 

21.62 

«• 

22.88 

25.54 

1 

22.90 

25.55 

3 

22.88 

25.53 

1 

10.15 

12.69 

15 
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Wavelength 

Element 

Multiplet 

(X) 

5041.063 

Si 

II 

Number 

5 

5056.020 

SI 

II 

5 

5105.541 

Cu 

I 

2 

5688.205 

Na 

I 

6 

5682.633 

Na 

I 

6 

5889.953 

Na 

I 

1 

5895.923 

Na 

I 

1 

5957.612 

Si 

II 

4 

5978.970 

Si 

II 

4 

6154.225 

Na 

I 

5 

6160.747 

Na 

I 

5 

6243.36 

A1 

II 

10 

6347.091 

Si 

II 

2 

6371.359 

Si 

II 

2 

6562.817 



1 


Excitation Potential 

Relative 

Low 

High 

Intensity 

(visual) 

10.02 

12.47 

15 

10.03 

12.47 

20 

1.38 

3.80 

1 

2.10 

4.27 


2.09 

4.27 

I 

0.00 

2.10 

18 

0.00 

2.09 

16 

10.02 

12.09 

4 

10.03 

12.09 

10 

2.09 

4.10 


2.10 

4.10 

J 

13.02 

15.00 

1 

8.09 

10.03 

24 

8.09 

10.02 

22 

10.15 

12.04 

25 
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APPENDIX 1 
StrcAk Photography 

The general arrangement of the streak camera la shown In Fig. I-l. 


O- - > 






M 


The object, 0, to be photographed Is 
reflected off the rotating mirror, M, 
on to the 35 ran film strip, F. 

The developed film strip pres 
entation is shown in Fig. 1 - 1 . 

The viilocity of the exposure on 
the film strip may be expressed as: 

Vj a 2nfR (I-l) 

where f a frequency of 

rotating mirror 

R a radius of curva¬ 
ture of the streak 
camera film holder. 

If the velocity of the luminous shock front, v, is desired at point 
A shown in Fig. 1-2, a tangent (i.e. slope or velocity) at that point may 
be constructed as shown in Fig. 1-3. 

The right triangle shown has sides 
v^ and mv where m is the magnifica¬ 
tion of the streak photograph and v 
is the desired velocity of the lum¬ 
inous shock front at point A. There¬ 
fore: 

tan 0 (1-2) 



1 

m '^f 



Combining 1-1 and 1-2, 

V . MB tan 9 (1-3) 

" Fig. 1-3 

R was determined by direct measurement on the streak camera. The 
quantities 8 and m were determined from the film strip, whereby f was obx« 
tained through the use of a counter and timer connected to the axis of the 
rotating mirror. 
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APPENDIX II 


Spectroscopic Techniques 


It was necessary to make a number of calibrations for the experimental 
equipment. Since the Gaertaer L 254 spectrograph and associated equipment 
will be available at the Postgraduate School for further research work, the 
results of the calibrations made are listed below for future reference. 

A. Silt width calibration 

The Jaws of the silt were checked and found to be parallel and In- 
excellent condlton. The various silt width settings were measured 
and plotted. Fig. II-l Is a curve of silt width setting versus 
the measured silt width. 

B. Timer Calibration 

The following timer settings were determined: 


Dial Setting 
(sec) 

0.5 

1.0 

1.6 

2.0 

4.0 


Measured Time 
(sec) 

0.5 
1.0 
* 1.5 
2.4 
4.9 


C. Step Filter Calibration 

The source of the transmission values for the step filter was 
unknown. Measurements made showed considerable variation with 
the rated transmission values provided. 


Step number, n 

Rated 

Measured 

1 

1.0 

1.0 

2 

0.627 

0.406 

3 

0.390 

0.244 

4 

0.250 

0.147 

5 

0.156 

0.065 

6 

0.100 

0.030 

7 

0.064 

0.012 
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D. Other Calibrations 




The spectrograph was cleaned, aligned, and focused by Aerospace 
Corporation after receipt of the Instrument from HOTS, China Lake, 
California. The Standard Lsmp and associated equipment was simi¬ 
larly calibrated by Aerospace Corporation prior to experimental 
measurements. 

The measurement of radiant energy spectroscopically Is a very difficult 
and sensitive problem. There are many factors Involved, which. If not pro¬ 
perly compensated for, will lead to meaningless values. In general, the 

2 3 A 

recommended methods of Sawyer ; Brode ; and Harrison, Lord and Loofbourow 

were followed In the Intensity determinations. 

Particular attention was given to the following Items: 

1. Pre-measurement checks 

2. Plate exposure 

3. Plate developing 

4. Plate calibration. 

Pre-measurement checks Included careful warming and Insertion of the 
spectroscopic plate, check of spectrograph alignment, focus and optical 
set-up, and the complete darkening of the room prior to measurement. 

Plate exposures of the shock tube firing and the Standard Lamp were 
taken through Identical optical set-ups with the exception of the Standard 
Lamp window which necesiltated a separate compensation In the data reduction 
calculations. The exposures were taken over a short time Interval to reduce 
the affect of ambient temperature and pressure variations on the spectrograph. 

Temperature and time criteria were closely adhered to in the developing 
and processing of the plates. The chemicals used were freshly made, and 
constant agitation was also employed. 

The densitometer used in plate calibrations was carefully checked to 
Insure that It was fully warmed and steadied out prior to use. After each 
densitometer run was made, a check was Immediately made to Insure that drift 
has not occurred during the run. 


AS 
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